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Guozhu Sun, 1 ’ 2 ’ 3 ’!^] Xueda Wen, 4 Ming Gong, 3,5 Dan-Wei Zhang, 5 Yang 
Yu, 5,2 Shi-Liang Zhu, 5 Jian Chen, 1 Peiheng Wu, 1 ’ 2 and Siyuan Han 3 

1 National Laboratory of Solid State Microstructures and Research Institute of Superconductor Electronics, 
School of Electronic Science and Engineering, Nanjing University, Nanjing 210093, China 
2 Synergetic Innovation Center of Quantum Information and Quantum Physics, 

University of Science and Technology of China, Hefei, Anhui 230026, China 
3 Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045, USA 
4 Department of Physics, University of Illinois at Urbana- Champaign, Urbana, IL 61801, USA 
5 National Laboratory of Solid State Micro structures, 

School of Physics, Nanjing University, Nanjing 210093, China 

Landau-Zener transition (LZT) has been explored in a variety of physical systems for coherent 
population transfer between different quantum states. In recent years, there have been various 
proposals for applying LZT to quantum information processing because when compared to the 
methods using ac pulse for coherent population transfer, protocols based on LZT are less sensitive to 
timing errors. However, the effect of finite range of qubit energy available to LZT based state control 
operations has not been thoroughly examined. In this work, we show that using the well-known 
Landau-Zener formula in the vicinity of an avoided energy-level crossing will cause considerable 
errors due to coherent oscillation of the transition probability in a single-passage LZT experiment. 

The data agree well with the numerical simulations which take the transient dynamics of LZT into 
account. These results not only provide a closer view on the issue of finite-time LZT but also shed 
light on its effects on the quantum state manipulation. 


Landau-Zener transition (LZT) has broad applications 
in atomic and molecular physics, quantum optics, con¬ 
densed matter physics, chemical physics, and quantum 
information science. For example, LZT has been applied 
to investigating the jump time and quantum Zeno and 
anti-Zeno effects of cold atoms in accelerated optical lat¬ 
tices mm, the behavior of molecular magnets at low 
temperature Pill], nonequilibrium phase transitions [5], 
and it is also exploited as a tunable beam splitter of wave 
functions to generate entangled multipartite states Eld- 
LZT also plays a key role in determining whether random 
optimization problems can be solved using the quantum 
adiabatic algorithm [8]. Recently, LZT’s potential for 
robust manipulation of coherent quantum states has at¬ 
tracted much attention in the context of quantum infor¬ 
mation processing EHE1I3HTSI because LZT may provide 
a simple and effective solution to the realization of high 
fidelity quantum state control without the need for pre¬ 
cise timing. 

The time-dependent Hamiltonian describing LZT in 
quantum two-level systems can be written in the generic 
form as 

H LZ {t) = -^e(t)a z -^Aa x , ( 1 ) 

where cr X:Z are Pauli matrices, e(t) = vt is the energy dif¬ 
ference between the two diabatic (crossing) basis states 
(i.e., the eigenstates | f) and I i) of the a z operator) 
controlled by an external parameter which depends lin¬ 
early on time t , and A is the constant gap between the 
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two instantaneous eigenenergy states |+) and |—) at the 
center of the avoided crossing e = 0, as depicted in 
Fig. 1(a). In such systems, when e{t) is swept through 
the avoided crossing, transitions between |±) with ener¬ 
gies E ± (t) = ±\^e(t) 2 + A 2 can occur and the transi¬ 
tion probability is given by the well-known Landau-Zener 
(LZ) formula 

Plz = , (2) 

where v = \de/dt\ is the Landau-Zener speed and we 
have set the reduced Planck constant K = 1. Equa¬ 
tion ([2| gives the probability of finding the system in 
the excited (ground) state at £/ = e(t —> +oo) when it is 
started in the ground (excited) state at e* = e(t —► —oo). 
By defining a = A 2 /Av as the adiabaticity parameter 
the LZ formula can be simplified to Plz = exp(—2-7ra). 
Although analytical solution to the problem cannot be 
obtained when £, and/or £f are finite, it is well known 
that for |ejj| A, the LZ formula provides an excellent 
approximation to the actual transition probability and 
Pi m Plz- However, when \eij\ A is not satisfied, 
the LZ formula may become quantitatively inaccurate or 
even qualitatively incorrect. In spite of some theoretical 
studies on the effects of finite \eij\/A on Pj,, there is an 
acute lack of adequate experimental evidence. 

On the other hand, understanding LZT with mod¬ 
erate values of ]^*,/|/A is in urgent need because this 
region of parameter space is important to quantum in¬ 
formation processing. For instance, in superconducting 
qubits the tuning range of energy level spacing is usu¬ 
ally limited to a couple of GHz or even as narrow as a 
few hundreds of MHz while A/27T could be as large as 
10 2 MHz [BU71CE9H26]. For quantum state control based 
on sweeping e through avoided crossings, understanding 
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LZT probability’s dependence on \eij |/A and the sweep¬ 
ing time is essential to high fidelity operation. The fi¬ 
delity of various techniques based on LZT relies critically 
on the accuracy of the LZ formula which predicts a sim¬ 
ple exponential dependence of Plz on the adiabaticity 
parameter a only. Therefore, for a qubit starting from 
the ground (excited) state the probability of finding it in 
the excited (ground) state after a single passage through 
the avoided crossing is assumed to be determined en¬ 
tirely by a (i.e., A 2 /4u) but not the detail of the process 
such as |e:jj|/A. Here, using an artificial atom — a su¬ 
perconducting phase qubit — coupled to a microscopic 
two-level system (TLS), we test the accuracy of the LZ 
formula in the region of |eq/|/A < 4.3. We show that in 
contrast to conventional wisdom, in the region of param¬ 
eter space most relevant to superconducting qubits, Pj, 
could deviate significantly from Plz determined by the 
LZ formula. Our experiment and numerical simulation 
demonstrate P± can oscillate coherently as a function of 
£f for constant a when \£ij\ is comparable to A, which 
is named as coherent Landau-Zener oscillation (LZO). 
Results 

In our experiment we use a superconducting phase 
qubit. However, since a single phase qubit does not have 
an intrinsic avoided energy-level crossing, we utilize an 
avoided level crossing arising from interaction between 
the qubit and a microscopic TLS [27] [28]. As discussed 
below in more detail, when the transition frequency of the 
qubit ui io is close to that of the TLS lutls , which is fixed, 
the first and second excited states of the coupled qubit- 
TLS system form an effective quantum two-level system 
described by the LZ Hamiltonian 0 . Note that to make 
quantitative comparisons between the theory/numerical 
simulation and the experiment without free parameters, 
all relevant system parameters, including the energy re¬ 
laxation and dephasing time of the qubit and the energy 
gap A, are obtained from direct measurements. 

A microscopic picture of the superconducting phase 
qubit is shown in Fig. 1(b). The qubit consists of 
an L « 770 pH superconducting loop intersected by a 
Al/AlOz/Al Josephson tunnel junction with a critical 
current I 0 « 1.4 /xA and a junction capacitance C ~ 240 
fF. By varying the magnetic flux applied to the supercon¬ 
ducting loop the potential energy of the qubit becomes 
asymmetrical. The ground state and the first excited 
state in the upper potential well, represented by |0) and 
|1) respectively, can be used as the computational basis 
states of the qubit. For an isolated qubit, the transi¬ 
tion frequency between |0) and 11), wio, is a single-valued 
function of the external flux bias <& x which is coupled in¬ 
ductively to the superconducting loop through an on-chip 
flux bias line. 

As shown in Fig. 2(a), however, the microwave spec¬ 
trum of the qubit wio(4 > a; ) has a rather large avoided 
energy-level crossing at « —2.8 m$o (with respect to 
the flux bias point at which u>iq/2tt ss 16.348 GHz) indi¬ 
cating significant interaction between the qubit and a mi¬ 
croscopic TLS [50 ! . The transition frequency between the 


TLS’ ground state |g) and excited state |e) and the qubit- 
TLS coupling strength are ujtls/ 2tt = 16.450 ± 0.002 
GHz and A/27T = 70.0 ± 0.5 MHz from the spectrum 
and vacuum Rabi oscillation, respectively. Note that in 
this coupled qubit-TLS system the time-dependent en¬ 
ergy difference between the two diabatic states involved 
in LZT is e(t) = coio(t) —u>tls which depends linearly on 
the flux bias to a good approximation. The relationship 
between the flux bias and e can be found from Fig. 2(a). 

Fig. 1(c) illustrates the experimental procedure used 
to observe coherent LZO. We begin by setting the ini¬ 
tial diabatic energy of the effective quantum two-level 
system £{ at about 100 MHz below ujtls with a static 
flux bias. The qubit is prepared in its ground state by 
waiting for much longer than the energy relaxation time 
Xj « 70 ns of the qubit. A microwave pulse is then 
applied to the qubit when it is biased at a fixed value 
£i/27r S3 —100 MHz. The microwave pulse coherently 
transfers the population of the qubit-TLS from |0c/) to 
one of the system’s eigenstates |—) through a process 
that is discussed in detail in Methods. The lack of os¬ 
cillation in Ti measurement taken at £; as shown in Fig. 
2(b) confirms that the initial state of the qubit-TLS sys¬ 
tem at t = 0 is indeed the eigenstate |—). As illustrated 
in Fig. 1(c), a time-dependent flux 4>(f) = §Lzt/t sp is 
then superimposed between t = 0 and t = t sp onto the 
static flux bias to sweep e linearly from ~ —100 MHz to 
its maximum value £/. The corresponding LZ speed v is 
thus ( £f—£i)/t sp . This is followed immediately by a 5-ns 
readout pulse which performs a projective measurement 
of the probability P± of finding the qubit in state |1) (i.e., 
the coupled system is in state |lg) corresponding to | j.) 
in Fig. 1(a)). 

We first measure P vs. t sp at a constant value of £f by 
keeping <&lz fixed while increasing t sp from almost 0 ns 
to 45 ns. The maximum t sp is selected to avoid too much 
influence of the qubit’s energy relaxation. By stepping 
$lz from 0 to —11 m$o the value of corresponding £f 
is then varied from about —1.4A to ~ 4A , in which the 
condition |e^| A is no longer satisfied. This procedure 
is repeated at each £f to obtain P±(£f,t sp ). Fig. 3(a) 
shows the dependence of P± on £/ and t sp . It can be 
seen that P± vs. t sp decays exponentially for £f < 0 
(d’iz € [—2.8 m$o, 0]) with a characteristic time Ti 
due to energy relaxation. As £/ becomes positive, P± 
vs. t sp becomes oscillatory. Since the avoided crossing is 
traversed only once, the observed oscillation in P± vs. t sp 
with constant £/ must be a consequence of the moderate 
value of £//A and is not caused by the Landau-Zener- 
Stuckelberg interference which requires multiple passages 
through the avoided crossing. 

It is worth noting that the observed oscillation is not 
a consequence of ill-prepared initial states with non- 
negligible probability amplitude in the excited state |+) 
of the effective Hamiltonian ([!]) because the microwave 
pulse used to initialize the system resonantly couples |0g) 
to |—) and has negligible coupling to |+) due to large 
frequency detuning. Furthermore, this process of trans- 


3 


ferring the system to the desired initial state |—) via a 
resonant microwave pulse is robust in the sense that it 
does not depend sensitively on the accuracy of the pulse 
duration t-MW- Deviation in pulse duration simply leaves 
some probability amplitude in 10 g) which has no effect on 
LZT other than reducing the visibility of the oscillation 
(see Methods for detail on the initial state preparation). 
Therefore, we are confident that the oscillation observed 
in the non-adiabatic region of the parameter space arises 
neither from Landau-Zener-Stiickelberg interference nor 
unwanted probability amplitude of |+) in the initial state. 
This is also supported by the good agreement between 
the results of experiment and numerical simulation shown 
in Fig. 3(b), which uses |—) as the initial state at the 
start of the single passage sweep. 

Discussion 

By replacing Ef with E t + vt, solving the problem of 
sweeping e in a finite range is transformed to finding P j, 
at finite time. Previous studies have discovered that LZ 
transition probability reaches the asymptotic value given 
by equation p]) at t t^z (assume e = 0 at t = 0), where 
tLz = 2- v /a/Amax(l, y/a) is called the Landau-Zener 
time, and oscillates in the vicinity of avoided crossing 
(corresponds to t < tLz) due to the transient dynamics 
mmm- Since corrections to the standard LZ formula 
are significant only if the adiabaticity parameter a < 1, 
the region of e within which transient dynamics plays an 
important role is given by £/ < tj p , where “<” means 
less than or comparable to. By examining the experi¬ 
mental data we find that the region of most noticeable 
coherent LZO coincides with £/ < t~ p which agrees well 
with the result of numerical simulation. These results un¬ 
ambiguously show that coherent LZO is originated from 
the transient dynamics of the LZT. 

Such coherent LZO has little effect on the adiabatic 
evolution. Because in the true adiabatic regime, by 
definition the system always stays in the instantaneous 
ground state and no LZT could occur. In order to find 
the region of approximate adiabatic evolution in our ex¬ 
periment, it is necessary to modify the definition of the 
adiabaticity parameter to a' = ^ UTLS ~^°'> + A _ F or the 
adiabatic theorem to hold a' 1 is required. As shown 
in Fig. 3(a), the white dashed lines represent a' = 10. 
It is clear that there is no coherent LZO observed in the 
region a' 3> 1. 

In a popular analogy to optics an avoided crossing acts 
as an effective beam splitter, with a transmission coeffi¬ 
cient corresponding to Plz in the LZ formula, for quan¬ 
tum wave functions. This beam splitter analogy has been 
applied successfully to the visualization and explanation 
of the behavior of superconducting and semiconductor 
qubits IS 0 H91 EDI G5H35]. In this analogy, a single 
sweep through the avoided crossing is equivalent to pass¬ 
ing a beam of light through the beam splitter only once. 
When |ejj| A, P± ss Plz and thus a greater LZ speed 
corresponds to a higher transmission coefficient of the 
beam splitter according to the LZ formula. But when 
\sij\ A is not satisfied, Pf differs greatly from Plz- 


As an example, P± vs. t sp , and thus the LZ speed v, with 
Ef/= 200 MHz is shown in Fig. 4(a). The maximum 
in the difference 5P between the experimental and 
those obtained from the LZ formula (2), shown in the 
inset of Fig. 4(a), can reach 0.21. The observation of co¬ 
herent LZO strongly suggests that when \£ij\ A is not 
met corrections to the LZ formula should be considered 
to avoid conceptual difficulties. 

Coherent LZO also has significant consequences on the 
coherent manipulation of quantum states of single qubits 
and coupled two-qubit systems based on LZT mm- 
For this approach of quantum state control, the LZ tran¬ 
sition probability Plz plays a central role since each sin¬ 
gle passage through the avoided crossing results in a uni¬ 
tary operation Ulz given by )12| 

Wl - P LZ e~ iVs -y/p^ \ 

VP^ —i\/l — P LZ e iV: ) ’ 

where ip s is the Stokes phase E3, which has no effect on 
the single-passage LZT process discussed here and thus 
can be set to zero for the sake of convenience. As men¬ 
tioned above, the transition frequency w 10 of most arti¬ 
ficial atoms, in particular the superconducting qubits, 
is limited to a couple of GHz. Because the speed of 
two-qubit operations is proportional to the inter-qubit 
coupling strength A, increasing \eij\/A by reducing A 
is undesirable. Hence, evaluating Ulz according to the 
LZ formula ([2]) could result in significant errors when 
Ef A is not satisfied. In order to conduct a quanti¬ 
tative analysis, the experimental data along the yellow 
dashed line corresponding to constant v. in the Fig. 3(a) 
are extracted and shown in Fig. 4(b). Oscillation in P|, is 
clearly observed and it is qualitatively different from the 
exponential decay predicted by equation ([2]), when deco¬ 
herence is taken into account. Suppose the initial quan¬ 
tum state is | 4-)- Then after a single passage through 
the avoided crossing, if one replaces Plz with Pf as in 
the asymptotic situation, the deviation SPf = Pf — Plz 
would be quite large. For example, when t sp = 12.8 ns 
the deviation SPf = 0.229, which is unacceptably large 
for coherent quantum state transformation. 

In conclusion, we have investigated the effect of finite 
energy (e) sweep (or equivalently finite time) on LZT 
probability Pj, experimentally. Single-passage technique 
is used to isolate the effect of finite Ef on Pj from that 
of interference caused by passing the avoided crossing 
multiple times. We find that P±(Ef/A,a = const) os¬ 
cillates when Ef is comparable to A and a < 1. The 
good agreement between the experiment and numerical 
calculation strongly supports the notion that coherent 
LZO is caused by the underlying transient dynamics of 
the finite time LZT which cannot be described by the LZ 
asymptotic formula. In this region of the LZT parameter 
space, corrections to the LZ formula must be taken into 
account, otherwise it will lead to substantial errors in 
quantum state operations based on LZT. The result also 
shows that when applying the simple beam splitter anal¬ 
ogy one should not automatically assume that greater a 
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(i.e., faster sweep) corresponds to larger transmission co¬ 
efficient (i.e., greater Plz) as implied by the asymptotic 
LZ formula. 

Methods 

Initial state preparation 

We first derive an analytical result explaining the lack 
of oscillation at the very beginning of LZT. The Hamil¬ 
tonian of the qubit-TLS system coupled to a microwave 
field is given as: (in the basis {| 0 g), |1 g), | 0 e), |le)}) 


H = 


0 fl m cos cot 

f l m COS Cot CO 10 

0 A/2 

0 0 


0 0 

A/2 0 

cotls flm cos cot 
rim cos U}t Wio + IOTLS 


(3) 

where fl m is the Rabi frequency, to is the microwave fre¬ 
quency, lotls is the energy difference between the ground 
state | < 7 ) and the excited state |e) of TLS, A is the cou¬ 
pling strength between the qubit and TLS, 6 = com — 10 
and S r = cjtls ~ wio are detunings. By rotating the 
frame, Hamiltonian ([3| can be transformed to the fol¬ 
lowing time-independent form [ 29| l30| 


Hi = 


1 

2 


- 5 r -25 n m 0 0 \ 

fl m —Sr A 0 I 

0 A S r fl m I 

0 0 fl m S r + 26 J 


(4) 


Next, we rewrite Hi in which the subspace spanned by 
{| 1 < 7 ), | 0 e)} is diagonalized: 


/ s r +26 _n _ a_q n \ 

/ 2 41V_ lira 4J¥ + u \ 

_^_f) _ 1 ,/A 2 4 . X 2 n _A_n 

4 N_* c ‘ rn 2 V a ' u r u 

_4_o n I./A 2 4- 6 2 _Wo 

4/V_|_ ^ Lrn u 2 V a rO r 4iV_ ^ Lrn 


(5) 


where N T = sj A 2 /4 + {-S r /2 ± v/z^TW/2) 2 - The 
basis for H 2 is now {|0g), |—), |+), |le)}. Note that in 
our experiment, before turning on the microwave the 
state is at ^(t = 0) = |0g). By turning on the mi¬ 
crowave {fl m ^ 0 ), | 0 <jr) is coupled to both |—) and 
|+). The resonance between |0 g) and |±) occurs when 
_Sr± 2 S _ _)_ A ' 2 + Sr/ 2 , from which we obtain the res¬ 

onant condition: 


— 1 LS ± ^ \/( u tls ~ ^io ) 2 + A 2 = A±. ( 6 ) 


In the limit of com ~ utls A, we have 10 = wio, 
which corresponds to the usual two-state Rabi oscilla¬ 
tion. Note that in this limit there is also a solution 
co = cotls- However, in this case the coupling strength 
is — 2 —> 0. The reason is that although the mi¬ 
crowave frequency could match that of TLS, coupling 
between the microwave and TLS is negligible which is 
confirmed by the absence of Rabi oscillation between the 
two states of the TLS in a separate experiment. In the 
other limit of com~ wtls = 0 , we have co = — + “ TLS ±A, 


and the dynamics have been thoroughly studied in Ref. 

m- 

In our experiment, we have (cotls ~ wio)/ 27 t ps 100 
MHz and A/27T ps 70 MHz, which means LZT occurs 
in the region where ( cotls ~ w io) ~ A. Because the fre¬ 
quency of the applied microwave is co = A_, which can be 
determined from the measured energy spectrum shown 
in Fig. 2(a), |0g) is resonantly coupled to |—), which 
is the eigenstate of H b = wi 0 |lff)(lff| + wrLs|0e)(0e| + 
Y(|l<?)(0e| + |0e)( 1 < 7 |). Although there is in principle also 
a coupling between | 0 g) and |+), the effective coupling is 
much smaller because of the large detuning, as discussed 
below. 

For A/27T ps 70 MHz, ( cotls ~ wio)/27t ps 100 MHz, 
the resonance between | 0 g) and |—) occurs at S/ 2 tt ~ 11 
MHz, we obtain N_/2tt ~ 36.7 MHz and N+/2tt ~ 116.4 
MHz. In our experiments, the coupling strength between 
| 0 < 7 ) and |—) is about 20 MHz and that between \Qg) and 
|+)would be x 20 MHz ~ 6.3 MHz. Because 6.3 MHz 
is comparable with 20 MHz, one may think that coupling 
between | 0 < 7 ) and |+) cannot be neglected. However, 
there is also a large detuning of about 122 MHz between 
10g) and |+). Therefore, the effective coupling between 
\0g) and |+) is reduced to (6.3 2 /122) ps 0.33 MHz and 
thus can be safely neglected. To be more precise, we cal¬ 
culated the population P± (where P± is the population of 
state |±)) after the application of a n pulse numerically, 
and it is found that P + /P_ ~ 5 x 10 5 . Based on this 
analysis, when co = A_, the dynamics can be described 
by the Hamiltonian in the subspace {|0g), |—)}: 


S r +2S 


H, = 


_A-O 

4 


_ A O _I,/A2 4- S 2 
4N- * Lrn 2 V a 1 u r 


(7) 


At the resonance — A±M = _ i ^/A 2 + S%, H 3 becomes 

where I is a 2 x 2 identity matrix. 

For initial state ^(f = 0) = |0g), the amplitude of |—) 
is 

C- = i sin -^-fl m t. (9) 


Considering |—) = — 5 ^ 7 1 Iff M— <? '~ + 2 jv^ + ~|0e), the am¬ 
plitude of 11 ) is 

A . A 

C\ — Ci g + C i e ~ Ci g — — i sin ^ fl m t ■ 

Therefore, with a microwave pulse of duration t^w 
which is used to prepare the initial state, we have 

P i = \Ci\ 2 ~ ^2 (sin ^-flmtMw^j ■ (10) 
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This is the reason why in the experiment we observe a 
usual Rabi oscillation instead of Rabi beating (30i which 
is indicated by the red circles, as shown in Fig.2(b). In 
addition, when microwave is turned off, the subspace 
{|+), |—)} is isolated from |0g) and |le). Projected into 
the subspace {|+),|—)}, the system stays in the eigen¬ 
state |—). This explains why we observe a monotone 
decay of P ( with no oscillations, as indicated by the blue 
triangles in Fig.2(b). 

Effect of tMW and £f/ A on LZT 

In this section, we discuss two factors that may af¬ 
fect the LZT probability, i.e., the width of the microwave 
pulse tMW used to prepare the initial state at and the 
end of the normalized diabatic energy sweeping £//A, 
respectively. 

After a microwave pulse, by projecting into the sub¬ 
space {|It/), |0e)}, the system is in the eigenstate |—). 
Then the dynamics of LZT can be described by Hb with 
a time-dependent wio(i), i.e., 


(11) 

To investigate the Landau-Zener diffraction effect, we 
sweep wio(t) across ojtls , he., 


wio(t) = uj 10 (t = 0) + iff, (0 < f < t sp ). (12) 


When wio(t = t sp ) — ojtls A, we expect that the 
Landau-Zener asymptotic formula holds and P± ss Plz- 
Thus no oscillations should occur in P±. This is confirmed 
by the result of numerical simulation shown in Fig. 5(c), 
where £f/2n = 1450 MHz « 20.7A, reproducing the ex¬ 
ponential decay behavior described by the asymptotic LZ 
formula independent of the initial state of the qubit-TLS 
system. In this case, the population in the qubit state 
|1) can be expressed as 

P i ^2 ( Shl X e-^ X e -WMw+i ap ) 

(13) 

where the first term reflects the effect of microwave dura¬ 
tion tMW iu preparing the initial state, the second term 
corresponds to the LZT probability, and the third term 
represents the relaxation effect. However, as 0 Jin(t = t sp ) 
moves towards ojtls , the situation uiio(t = t sp )— ojtls 
A does not hold any more, and we observe oscillation fea¬ 
tures in the t sp direction, as shown in Fig. 5(a) (experi¬ 
ment) and 5(b) (numerical simulation). Notice that Fig. 
5(a) and 5(b) also confirm that the effect of imprecise 
7 r pulse is an incomplete transfer of system from |0g) to 
|—), which reduces the probability amplitude of |—) from 
the maximum value, instead of resulting in non-negligible 
probability amplitude in the unwanted |+) state. 
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FIG. 1: Circuit and experimental procedure, (a) A gen¬ 
eral avoided energy-level crossing with diabatic basis states 
(the dashed lines) and adiabatic basis states (the solid lines). 
The constant gap A between the two instantaneous eigenen- 
ergy states at e = 0 is the tunneling amplitude, i.e., the cou¬ 
pling strength, (b) Optical micrograph of the sample with 
Al/AlO^/Al Josephson junctions on the silicon substrate. 
MW, and RO represent the microwave pulse to pre¬ 

pare the initial state, the sweeping flux bias to induce LZT 
and the 5-ns flux bias pulse to tilt the potential well in or¬ 
der to readout the qubit state, respectively. The planar coils 
and dc-SQUID magnetometer are coupled inductively to the 
qubit, (c) A typical time profile of the manipulation and mea¬ 
surement waveforms employed to perform the single-passage 
Landau-Zener experiment. 


a S! 2*= (co.^mra) / 2ji (MHz) b 
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FIG. 2: Spectroscopy and Rabi oscillation, (a) Mi¬ 
crowave spectroscopy measurement of the coupled qubit-TLS 
system. The splitting is A/2n = 70.0 ± 0.5 MHz centered 
at uj/2n = 16.450 ± 0.002 GHz. The beginning point of the 
flux bias for the single passage LZ sweeping is denoted as 0 
m<I>o, corresponding to £i/2i: = (gjio — ujtls)/2.-k « —100 
MHz in the upper abscissa, (b) Rabi oscillation and Ti at Ei, 
respectively. The experimental data (the red circles and blue 
triangles) agree well with the theoretical fits (the black solid 
lines). 
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FIG. 3: Coherent LZO. (a) Experimentally measured Pj_ 
vs. £/ and t sp . (b) Numerically calculated Pj, vs. e/ and 
t sp with all input parameters obtained from the experiment. 
The white dashed lines correspond to the value of modified 
adiabaticity parameter a' = 10. Notice that in the region 
below (above) the lower (upper) white dashed line, one has 
a! >> 1 thus the system evolves adiabatically and no oscilla¬ 
tion in Pj_ is expected as confirmed experimentally. The LZ 
speed v equals the slope of any straight lines originated from 
the lower-left corner of the t sp -£f plane. For example, the yel¬ 
low dashed line in (a) has v = 400/19.5 « 20.5 MHz/ns. For 
the sake of clarity, the temporal evolution of the system along 
the yellow dashed line (constant v) is presented separately in 
Fig. 4(b). 
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FIG. 4: P_i oscillation at constant LZ speed, (a) Mea¬ 
sured Pi as a function of t sp (the red squares) with Sf/2n = 
200 MHz which clearly shows oscillation in the region of a < 1 
which compares well with the result of numerical calculation 
(the solid line). This is in stark contrast to the smooth ex¬ 
ponential decay expected from the LZ formula (the dashed 
line). The inset is the difference SP\ = Pj_ — Plz- (b) The 
measured (the red squares) and numerically calculated (the 
solid line) P; vs. t sp with constant LZ speed v « 20.5 MHz/ns 
corresponding to evolving along the yellow dashed line in Fig. 
3(a). Again, Pj. oscillates in the region where the adiabatic 
condition a > 1 is not satisfied which is not expected from 
the LZ formula (the dashed line). 
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FIG. 5: The effect of tMW on finite time LZT (a) Ex¬ 
perimentally measured and (b) numerically calculated Pj, vs. 
tMW and tap for £f/2n = 200 MHz showing that the effect of 
imprecise n pulse is to reduce the visibility of the oscillation 
of Pj. vs. tap by reducing the probability amplitude of the 
desired |—) state. Because the microwave pulse is resonant 
with the |0g} |—) transition while largely detuned from 

the | Op) |+) transition even a significant deviation from 

a 7r pulse would only result in negligible transfer of popula¬ 
tion to |+). Furthermore, since |0<?) does not participate in 
the single passage LZ process, the observed oscillation could 
neither be due to LZS interference nor non-negligible popula¬ 
tion in |+} at the beginning of each e sweep. For comparison, 
we also present the numerically calculated P\.(tMW,t sp ) for 
£//27 t = 1450 fa 20.7A in (c). The result shows the exponen¬ 
tial decay behavior described by the asymptotic LZ formula 
as expected for £//A 3> 1. 


